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i?e(ylo), i?e(A2) and RG evolution for Nf = 3* 

Keunsu Choi*^ and Weonjong Lee*^ 

^School of Physics, Seoul National University, Seoul, 151-747, South Korea 

We present results of Re{Ao) and Re{A2) calculated using HYP staggered fermions on the lattice of 16'^ x 64 
at /3 = 6.0. These results are obtained using leading order chiral perturbation in quenched QCD. Buras's original 
RG evolution matrix develops a removable singularity for Nf — 3. This subtlety is resolved by finding a finite 
solution to RG equation and the results are presented. 



In the standard model, we can express Re{Ai) 
in terms of the matrix elements of the effective 
weak Hamiltonian between hadronic states as fol- 
lows: 

Re{Ai) - 9£\v^^Vus\[j2 z,{fi){Q,{fi))i 
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-fi?e(T)^y,(^)(Q,(^))j 



(1) 



1=3 



where Vij is an element of CKM matrix. Here, 
Zi{^) and UiifJ.) are the Wilson coefhcients which 
are obtained analytically. The hadronic matrix 
elements {Qi{n))i are defined as 



{Q^{^l))I = {nni\Q,{ti)\K) 



(2) 



Since these matrix elements are in a highly non- 
perturbative regime of QCD, it is necessary to 
calculate them using non-perturbative tools such 
as a numerical method based on lattice gauge the- 
ory. In fact, we calculate {n\Qi\K) and {Q\Qi\K) 
on the lattice using HYP staggered fermions and 
we construct {ttt: i\Qi\K) out of them using chiral 
perturbation in its leading order p. The Re{T) 
in Eq. |^ is defined as 



i?e(r) 
A/ 



^i?e(At/A„) = 0.002 



= VfsVfd 



In Eq. (P), Qi (i = 3, 4, • • • , 10) represents the 
QCD and electroweak penguin operators. Since 

* Presented by W. Lee. Research supported in part by 
BK21, by Interdisciplinary Research Grant of Seoul Na- 
tional University, and by KOSEF contract ROl-2003-000- 
10229-0. 



ZijUi ~ 1, the contribution from the penguin op- 
erators to Re{Aj) is suppressed by Re{T). There- 
fore, in practice, only the current-current opera- 
tors {Qi) {i = 1,2) can contribute dominantly 
to Re{Ai) whereas the contribution from pen- 
guin operators is negligible. In the case of the 
direct CP violation, only penguin operators can 
contribute to e'/e and current-current opera- 
tors do not play any role in it. In this regards, 
we can say that e'/e and Re(Ai) are completely 
different physics. 

In Fig. ^ we show individual channels con- 
tributing to Re^Ao) with the total sum in the 
last column. The dominant contribution comes 
from (Qi) {i = 1,2) as expected. Each chan- 
nel is obtained by fitting the data to a function: 
/i(mx) — Co + cim^ and extrapolating to the 
chiral limit. We use the matching formula given 
in 3 4 to convert lattice results into continuum 
values in the NDR scheme at the scale of /i — 1/a. 
By assuming that these values are close enough 
to those for Nf = 3, we run them from fi — l/a 
down to /i = TOc, using the RG evolution equation 
for Nf = 3. 

There are two alternative methods to tran- 
scribe the Qe and Q5 operators to the lattice: 
the standard (STD) method and the Golterman 
& Pallante (CP) method 0. In the case of e'/e, 
the dominant contribution to the A/ =1/2 am- 
plitude comes from (Qg) and so there is a big dif- 
ference in e'/e between the STD and GP methods 
12 ■ In the case of i?e(Ao), however, the contribu- 
tion from (Qe) is extremely suppressed by i?e(r) 
and so both methods do not make much difference 
to Re{Ao). 
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Figure 1. Re{Ao), linear fitting for chiral extrap- 
olation, combined at /i = TOc. 



Figure 2. Re{Ao), linear fitting for chiral extrap- 
olation, combined at fi — 1/a. 



In Fig. 13 we show the same kind of plot as in 
Fig. ^except for the scale at which the Wilson co- 
efhcient and the matrix elements are combined. 
We set the scale to /i = 1/a in Fig. |21 and to 
jjL = rric (charm quark mass) in Fig. ^ The con- 
tribution from each individual channel can fluc- 
tuate depending on this matching scale but the 
total sum should be invariant. 

In Fig. O we show the same kind of plot as 
in Fig. n except for the fitting function. Each 
channel in Fig. O is obtained by fitting the data 
to a function: /2(m-K) = cq -I- cim\ + C2mj^ and 
extrapolating it to the chiral limit. Comparing 
Fig. n] and |21 we observe that the dependence of 
the chiral extrapolation on the fitting function 
is large. These results are consistent with the 
experimental value of i?e(Ao) = 33.3 x 10~^ GeV 
within the systematic and statistical uncertainty. 

In Fig. 21 we present individual channels con- 
tributing to Re{A2) and the total sum in the last 
column. Each channel is obtained by fitting the 
data to a linear function /i(TOi^) and extrapolat- 
ing it to the chiral limit. The result of Fig. 0] 
is slightly lower than the experimental value of 
Re{A2) = 1.5 X 10-8 GeV. 

Now, let us address a subtle issue on the RG 
evolution for Nf = 3. The AS = 1 effective 
Hamiltonian can be written as follows: 



n 



^ 10 
^ i—1 



(3) 



The RG (renormalization group) equation for the 
Wilson coefficient Ci is 



(4) 



where (3{g) is the QCD beta function: 



Pig) = 



9^ 



Pi 



9" 



Pie ,,^_9X9 (5) 



■I67r2 '^'(167r2)2 ^""(167r2)2 
and 7(5, a) is the anomalous dimension matrix: 
7(5, «) - 7s(.9') + ^r(g2) + ... (6) 



+ 



(7) 



^'^y > Att ^' (47r)2 

A solution to the RG equation can be expressed 
in terms of the evolution matrix. 

C(A*) = C/(/i,/i',a)^(/i') (8) 

Here the RG evolution matrix is, in general. 



(9) 



Iff ^ T ( r'"^'. 

\Jg(ni2) Pig) 

In the perturbative expansion, we can express the 
RG evolution matrix as follows: 

U{mi,m2,a) — U{mi,m2) + —R{mi,m2) (10) 

47r 

where C/(mi, 1112) represents the pure QCD evolu- 
tion and i?(mi,TO2) describes the additional evo- 
lution in the presence of the electromagnetic in- 
teraction. The pure QCD evolution matrix which 
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Figure 3. -Re(Ao), quadratic fitting for chiral ex- 
trapolation, combined at fj, — rric. 



Figure 4. Re{A2), linear fitting for chiral extrap- 
olation, combined at n = m^. 



Buras, et al. provided originally |6I7| . is 

as(mi) 



[/(mi, 7712) 



1 



4tt 

as{m2) 
47r 



J [/(°)(mi,m2) 



(11) 



where [/('^^(toi, TO2) denotes the evolution in the 
leading logarithmic approximation and J corre- 
sponds to the next-to-leading correction. 
The J matrix can be expressed as follows: 

= VSV-^ , G = V-^-i'^P^V (12) 
(0) Pi 



J 



Sij — 



2/3o + if^ - 7 



(13) 



where 1^ is a matrix which diagonalize the 7^ 
matrix and 7!°^ represents corresponding diago- 
nal elements. Here, note that the J or 5 ma- 



Jo) 



7 



(0) 



0. In fact, 



„(o) 



7j ' =0 when i — 



trix is singular if 2/3o + 7, 
for Nf = 3, 2/3o + 7^ ^ 

j — 7. However, since U{mi,m2) is not singular, 
this singularity must be removable at the next- 
to-leading order. After we make a correct combi- 
nation at the next leading order and remove the 
singularity, we can make the U {mi, 1112) finite. 

C/(mi,m2) = UQ{mi,m2) + ^VA{mi,m2)V^^(U) 

At: 

where, if aj ^ ai + I (note that Oj = 7j^"-'/(2/3o)), 



A{mi,m2) = Si. 



a. 



,{m2)Y 
,(mi) J 



-as{m2) 



and if aj — at 
A{mi,m2) = 



as{m2) 
a.s{mi) 

1, 

-fas{m2) — 

, Q!s(to2) 

• m 



(15) 



I ^ / (16) 

Here, note that the A matrix is finite while J 
matrix is singular. 

In the case of the R{mi,m2), the same kind 
of removable singularity makes the RG evolution 
much more complicated, which we will present in 
Ref. [S|. 
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